%dips bepalen ten gevolge van een tijdelijke lijnfout

clear all
tic
einde_interval = 1;   %# lijnen bekijken
teller = 0;
alle_spanningen = zeros(200*einde_interval,1);
tijd = zeros(200,einde_interval);
for k=1:einde_interval;
    
    %clear all;
    
    filename = sprintf('bus%03d',k);
    initpsat;
    Settings.plottype = 2;
    Settings.freq = 60;
    Settings.absvalues= 'on';
    runpsat(filename,'pf');
    Settings.fixt = 1;
    Settings.tstep = 1;%0.05; %TRAAAAAAAAAAAAAG
    Settings.plottype = 2;
    Settings.freq = 60;
    Settings.absvalues= 'on';
    
    filename = sprintf('bus%03d.m',k);
    
    idx_begin = DAE.n+Bus.n+1;
    idx_eind = DAE.n+2*Bus.n;
    
        fid=fopen(filename,'a');
        
        Varname.idx(:,1) = 0;
        teller1 = 1;
        for i = idx_begin:idx_eind
            Varname.idx(teller1,1) = i;
            teller1 = teller1+1;
        end
        
        fclose(fid);
    
    
    
%     runpsat('td');
%     %begint vanaf index 0... (STOM!!!)
%    v = genvarname('voltages_fout_lijn', who);
%    eval([v ' = Varout.vars']);%(:,(size(Bus.con,1)+1):2*size(Bus.con,1))']);
%    lengte = size(Varout.vars,1);
%    
%    %spanning per tijdsstap per bus
%    alle_spanningen(1+200*teller:200*teller+lengte,:) = Varout.vars;
%    
%    t = genvarname('tijd_fout_lijn', who);
%    eval([t ' = Varout.t']);%(:,(size(Bus.con,1)+1):2*size(Bus.con,1))']);
%    
%    %tijd per stap per lijn
%    tijd(1:lengte,k) = Varout.t;
   %closepsat;
   teller = teller+1;
   
%    %save dips tijd alle_spanningen;
%    save lijnfout tijd alle_spanningen einde_interval;
end

 %save dips tijd alle_spanningen;
   save lijnfout tijd alle_spanningen einde_interval;
% initpsat;
% Settings.plottype = 2;
% Settings.freq = 60;
% Settings.absvalues= 'on';
% runpsat('IEEE_14_bus_original','pf');
% Settings.fixt = 1;
% Settings.tstep = 1;%0.05;
% Settings.freq = 60;
%     Settings.absvalues= 'on';
% runpsat('td');
% normaal = Varout.vars;
toc